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Abstract 

The t — b — t unification with positive Higgs mass parameter /i in the minimal supersymmetric 

43 ! standard model prefers "just so" Higgs splitting and a light gluino < 500 GeV which appears to be 

Ph. 

ruled out by the recent LHC searches. We reanalyze constraints on soft supersymmetry breaking 

Q_|. 

. parameters in this scenario allowing independent splittings among squarks and Higgs doublets at 

(— ! . 

the grand unification scale and show that it is possible to obtain t — b — r unification and satisfy 
experimental constraints on gluino mass without raising supersymmetry breaking scale to very 
high value ~ 20 TeV. We discuss the origin of independent squark and Higgs splittings in realistic 
SO (10) models. Just so Higgs splitting can be induced without significantly affecting the t — b — T 
' unification in SO(10) models containing Higgs fields transforming as 10 + 126 + 126 + 210. This 

\q \ splitting arises in the presence of non-universal boundary conditions from mixing between 10 and 



o 

(N 



X 



other Higgs fields. Similarly, if additional matter fields are introduced then their mixing with the 
matter multiplet 16 is shown to generate the squark splitting required to raise the gluino mass 
within the t — b — r unified models with positive \x. 
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I. INTRODUCTION 



Grand unified theories (GUTs) based on SO (10) group not only unify the gauge inter- 
actions but also lead to a unified framework for the Yukawa couplings and hence fermion 
masses. In particular, SO (10) model with a 10-plet of higgs coupling dominantly to the 
third generation implies an equality yt = Ub = Ut of the t — b — r Yukawa couplings at the 
GUT scale. Quite independently, the renormalization group (RG) running of the Yukawa 
couplings in a softly broken minimal supersymmetric standard model (MSSM) can lead 
to the t — b — t unification at the GUT scale making the supersymmetric SO (10) broken to 
the MSSM at the GUT scale an attractive theory of unification. 

t—b—r unification at the GUT scale is however not the most generic property of the MSSM 
but follows only for a restricted set of boundary conditions for the soft supersymmetric 
breaking terms. These restrictions mainly arise due to the need of significant threshold 
corrections 0, to the b quark mass required for the t — b — r unification and difficulties in 
achieving the radiative electroweak symmetry breaking (REWSB) in the presence of large b 
and r Yukawa couplings js[ 0] . Both of these depend on the soft breaking sector. It is realized 
(5), @j that t — b — r unification generally requires departure from the universal boundary 
conditions assumed within the minimal supergravity (mSUGRA) framework. Universality of 
the gaugino masses is enforced by the SO (10) invariance if it is assumed that supersymmetry 
(SUSY) is not broken at the GUT scale by a non-trivial representation contained in the 
symmetric product of two adjoints of SO (10). In contrast, the soft masses mi 6 , m w for 
sfermions belonging to 16a/ and the Higgs scalars belonging to the 10h representations are 
allowed to be different and are also required to be so to obtain t — b — r unification. In 
addition to this SO (10) preserving non-universality, one also needs to introduce explicit 
50(10) breaking non-universality. Such non-universality can be induced spontaneously by 
a non-zero D-term (DT) which introduces splitting within 16a/ of squarks and 10^ of the 
Higgs simultaneously p, 0] . In several situations, one also needs to assume that only the 
MSSM Higgs fields H u , H d split at the GUT scale. This is termed as "just so" Higgs splitting 
(HS) @|. 

Restrictions placed on soft parameters by the t — b — r unification, the LEP and LHC 
bounds on the masses of the SUSY particles and other flavor violating observables have been 
worked out in detail in number of papers j8-13|. Two viable scenarios and their properties 



have been identified (see 14] for details and references therein). These depend on the sign of 



the n parameter of the MSSM. For example, one can achieve an exact t — b — r unification in 
mSUGRA itself for negative [i. But this needs very heavy SUSY spectrum with m ~ 5 — 12 



TeV and m\/2 ~ (1.5 — 2)mo [10] . Also, perfect t — b — r unification with relatively light 
SUSY spectrum (~ 2 TeV) can be obtained with the introduction of DT or purely Higgs 
splitting. This appears to be the best and testable scenario as far as the t — b—r unification is 
concerned. But the supersymmetric contribution to the muon (g — 2) is negative in this case. 



This adds to the existing discrepancy between theory 15j and experiments 16|. Scenario 



with positive \i proves better and allows the theoretical prediction for (g — 2) to agree with 
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experiments within 3a. In this case, achieving t—b—r unification becomes considerably more 



difficult. The mSUGRA in this case at best allows unification at 65% level [10|. Even with 
non-universal boundary conditions, one needs specific relations between the soft parameters, 
mio ~ 1.2mi6, Aq ~ — 2miQ and myi <C m\Q together with tan/3 ~ 50 in order to achieve 
t — b — r unification j^]. The DT splitting in this case, allows unification at most 90% level 
but this requires m\§ > 10 TeV and a gluino mass < 500 GeV. Just so HS works much 
better than DT splitting and leads to the perfect Yukawa unification for m 16 ~ 10 TeV but 
gluino is still light. In both these scenarios, the sparticle mass spectrum is characterized by 
lighter gluino which is within the reach of current LHC searches at \fs = 7 TeV, whereas all 
scalar sparticles have masses beyond the TeV scale. The light gluino mainly decays through 
a three body channel g — > bbx° leading to multijets plus missing transverse energy. The final 
states may also contain dileptons if t is lighter than b. Recently, the ATLAS experiment 
with 2.05 fb _1 data collected at y/s = 7 TeV has excluded the light gluino masses below 620 
GeV in 5O(10)+HS model As we will show later, this experimental limit on gluino 
mass rules out t — b — r unification better than 90% for m M ~ 10 TeV. The strong bound 
on the gluino mass follows from the need of suitable threshold corrections in bottom quark 
mass to achieve t — b — r unification which requires the hierarchy muz <C rnie. One can thus 
raise the value of m^j and hence the bound on the gluino mass by raising m^. The gluino 
mass can be pushed in this case beyond the present experimental limit but at the cost of 
choosing m 16 > 20 TeV 
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The best viable scenario with fi > corresponds to just so Higgs splitting and large 
m 16 and m 10 . Theoretically, both these features are unsatisfactory. A large SUSY scale is 
unnatural and just so HS breaks 50(10) explicitly. Just so HS can be indirectly introduced 
through the right handed neutrinos at the intermediate scale. Apart from causing problem 
with the gauge coupling unification, this case also does not do as well as the arbitrarily 



introduced just so HS, see [12j. We wish to discuss here possible ways to improve on both 
these aspects. Specifically, we show that just so HS arises naturally in realistic 50(10) 
models containing additional Higgs fields, e.g. the one transforming as 126, 126, and 210 
representations under 5O(10). Realistic fermion masses can be obtained if these fields are 
present and the SU(2)l doublets contained in them mix with each other. Moreover, all the 
angles involved in such mixing need not be small. We show that significant Higgs doublet 
mixing can generate just so HS in the presence of non-universal but 5O(10) preserving 
boundary conditions at the GUT scale without inducing any splitting between squarks or 
without significantly upsetting the t — b — r Yukawa unification. 

5O(10) models also allow an interesting possibility of matter fields mixing among them- 
selves [19|. This leads to just so squark/slepton splitting similar to the just so HS occurring 
due to Higgs mixing. Impact of such mixing was used earlier 19[ to obtain departure from 
the b-T unification that follows in SU(5) or 5O(10) models. Here we discuss phenomeno- 
logical implications of such mixing in the context of the t — b — r unification. We discuss 
explicit example leading to independent squark and Higgs splitting and show that the pres- 
ence of such splittings helps in raising the gluino mass prediction without raising the SUSY 



3 



parameters m\% to high values around 20 TeV. 

This paper is organized as follows. In the next section, we present a short review of the 
basic features of the t — b — r unification and discuss the existing phenomenological results. 
We also update the existing results incorporating the recent limit on B s — > from LHCb 



20|| . The viability oi t — b — r unified solutions in the presence of independent Higgs splitting 



and squark splitting (SS) are discussed in Section HTT1 In Section HVt we discuss how HS and 
SS can arise in the realistic versions of S'O(IO) models. The study is summarized in the last 
section. 



II. t-b-r UNIFICATION IN MSSM 



In this section, we review numerical and analytic results presented in the literature js 
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in the context of the t — b — r unification. We have re derived several existing numerical 
results in a way which optimizes the rate of the B s — > fi + n~ to make it consistent with the 



recent more stringent experimental [20| bound without pushing the SUSY scale to a higher 
value. Before discussing this, we summarize aspects of the t — b — r unification which allows 
understanding of the salient key features. 

The hypothesis of the t — b — r unification assumes that at the GUT scale 

m b m T m t 

Vb = a =Vt = ^ =Vt= 7—5 , (1) 

v cos p v cos (3 vsmp 

where v ~ 174 GeV. This equation is motivated by a simple 50(10) model containing only 
a single 10-plet Higgs. The appropriate choice of the free parameter tan f3 in Eq. (ED) can 
always allow equality of y t with y b or y T at the GUT scale. However it is well known [13[ that 
yb and y T and hence all three of them derived from Eq. ([1]) using the experimental values 
for fermion masses extrapolated to the GUT scale do not unify for any value of tan /3. The 
degree of unification of three Yukawas or lack of it is usually measured by the parameter 

_ Max.(y t ,y b ,y T ) 

Ktbr = ri — — — . (2 

Mm.{y t ,y b ,y T ) 

The parameter R t b T is defined at the GUT scale. Variation of R tbT with tan (3 obtained 
by using the tree level Yukawa couplings and fermion masses at Mz is shown in Fig. (CQ). 
For comparison we also show similar ratios R bT and R tr defined using only two of the 
Yukawa couplings. The extrapolation from Mz to the GUT scale is done using the 1-loop 
RG equations which are independent of the details of the soft SUSY breaking. Fig. ([Q) 
explicitly shows that the three Yukawas do not meet for any tan (3 and the reason is that b 
and r Yukawa couplings never meet. The best value of Rtb T seen from Fig. (CQ) is around 
1.2. Thus any scheme which tries to achieve t — b — r unification should do better than this 
tree level value. 

It is known that the tree level Yukawa couplings, particularly that of the b quark receive 0, 
0] significantly large radiative corrections once the supersymmetry is broken. The corrected 
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FIG. 1. Tree level Yukawa unification as a function of tan/?. Rij defined in Eq. ([2]) measures the 
closeness of j/j and yj at the GUT scale. 



y b can be written as 



y b = yt™ cos P(l + Ay 9 b + Ay- 



(3) 



where (...) contains the electroweak suppressed SUSY corrections, the standard model (SM) 
electroweak corrections and logarithmic corrections which are sub dominant. The dominant 
correction Ay b (Ay b ) induced by the gluino (chargino) exchange is approximately given by 



Ayl rs -^-/itan/3 



Ay2 



3tt 



y± 

167T 



/itan/3 — \- 
mf 
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(4) 



(5) 



where m 



b 2 



l m t 2 , 



is mass of the heaviest sbottom (stop). 



The presence of tan/3 makes the radiative corrections significant. The corrections to top 
Yukawa is inversely proportional to tan /3 while corrections to tau Yukawa is proportional to 
tan/3 but electroweak suppressed. In order to achieve unification, one needs to reduce R bT 
and hence y b compared to its tree level value by about 10-20%. This requires that Ay b + Ay b 
in Eq. ([3]) should be negative. Since the gluino induced contribution dominates over most of 
the parameter space, one can make the radiative corrections negative by choosing a negative 
fi. As a result, models with negative \x achieve t — b — r unification more easily. For positive 
/i, the chargino contribution has to dominate over gluino and it should be negative. This 
can be satisfied with a negative Aq and light gluino with \Aq\, mie ^> mi/2- As a result, all 
the scenarios of t — b — r unification with positive /i lead to a light gluino and very heavy 
SUSY spectrum as borne out by the detailed numerical analysis [8l-ll2 . 

The requirement of a light gluino as argued above directly conflicts with the requirement 
of the REWSB unless an explicit HS is introduced. This can be seen as follows. In large 
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tan/3 limit, the REWSB can be achieved if 



m 2 H > Mf/2 



Am 2 H 



m 



m% v > M z , 



(6) 



where d are soft scalar masses evaluated at the weak scale. Starting with a positive value 
at Mqut, m2 H u gets driven to large negative values by a large y t and the first equation gets 
satisfied. But the large ?/&, y T as required in the t — b — r unification drives even more 
negative and conflicts with the second requirement. In addition to the Yukawas, the gaugino 
and scalar mass terms also contribute to the Am 2 H . The former contribution is positive while 
the latter is negative, see the semi analytic solution of the 1-loop RG equations presented 
for example in 0]. In Fig. (j2J), we show the running Am^ for different values of m 1( / 2 . 
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FIG. 2. The solution of 1-loop RGE equation for Am], = m 2 Hd — m 2 Hu in mSUGRA for mo(= 
mio = m.16) = 1 TeV, Aq = and tan/3=50. 



As can be seen from the figure, second of Eq. (jSJ) can be satisfied by choosing rrii/2 > rriQ 
such that the gaugino induced contribution in Am 2 H dominates. On the other hand, small 
mi/2 and hence light gluino around m~ g < 500 GeV is required if significant corrections to 
is to be obtained in case of \i > 0. This corresponds to < 200 GeV for which REWSB 
cannot be achieved unless one introduces splitting between Higgs fields at the GUT scale 
itself. This is clearly seen from Fig. (j2J). Moreover, the case in which only Higgs splitting is 
considered is more favorable than the D-term splitting. This follows 0] qualitatively from 
Eqs. (jU EJ) which implies 



\M1 



llTT- 



m?, m 



A f \ 771- 



(7) 



One finds m 4 ~ 2 ~ for mSUGRA as well for just so HS. The D-term splitting introduces 
sfermion splitting together with Higgs splitting and leads to m t ~ 2 > m~ b2 . As a result, one 
needs to choose even a lighter gluino or a larger \A t \ to suppress the gluino induced correc- 
tions in yb. As we will show later in this paper, the additional squark splitting can instead 
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reduce the ratio m^/m^, and make it less than one. This allows significantly higher gluino 



mass. 



The above qualitative features are borne out by several numerical studies presented in a 
number of papers js-12]. A list of different scenarios proposed to achieve t — b — r unification 
for fi > is given in Table (1) in Ref. 12J. Among all the proposals, the 5*0(10) model 
with just so HS is found as the best scenario which leads to an exact t — b — r unification 
corresponding to R — 1. It is shown 
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that HS works particularly well for large m 16 
and the Yukawa unification R < 1.02 can be achieved if mi 6 > 10 TeV. We update this 
analysis for the following reasons. In lOj, the t — b — r unified solutions were obtained 
considering the experimental constraint BF(i? s — > yU + /i~)( cxp ) < 2.6 x 10 -6 . The recent data 
collected by LHCb 20[ experiment has improved this bound significantly. The current limit 
BF(B S -> /i+/i~) (cxp) < 4.5 x 10" 9 isjhree order of magnitude stronger than old bound. As a 



result, all the solutions obtained in 10[ are found inconsistent with new limit on B s — > fi + fi 



see, Tabel (1) in [lOj). We repeat the old analysis considering the new limits on B f — > 



and b — > 57. In addition, we also consider the present constraints on B — > rv T 
was not considered in the old analysis. 
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We use the ISASUGRA subroutine of ISAJET 7.82 [22j in our numerical analysis. For 
given boundary conditions (soft SUSY parameters at the GUT scale and m t , tan/3 at the 
weak scale), ISASUGRA solves full 2- loop MSSM RG equations and incorporates 1-loop 
SUSY threshold corrections in all the MSSM sparticles and in the masses of third generation 
fermions. Moreover, it checks for (a) non-techyonic solutions and (b) consistent REWSB 
using the minimization of one-loop corrected effective MSSM Higgs potential. Once these 
conditions are satisfied, we calculate R using Eq. fl2]). Then using CERN's subroutine 
MINUIT, we minimize R. Finally, we calculate branching factor for b —> 57, B s ^ A i+ /i~ 

For 



using the IsaTools package 
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and B — > rv T using the expressions given in 
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completeness, we also estimate the relic abundance of neutralino dark matter QcDMh 2 and 
the SUSY contribution to anomalous magnetic moment of muon Aa M (where a M = (g — 2)/2) 
using IsaTools. On the acquired solutions, we apply the following constraints obtained from 
experimental data: 



BF(B S -> < 4.5 x 10~ 9 

2.78 x 10~ 4 < BF(6 sj) < 4.32 x 10~ 4 (3a) 
0.62 x 10" 4 < BF(B ru T ) < 2.66 x 10~ 4 (3a) 
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21 
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(8) 
(9) 
(10) 



Further, we impose the mass bounds given in PDG 25[ on all sparticles including the LEP 
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bound on the mass of lightest Higgs (m^ > 114.4 GeV). We use m t = 172.9 GeV in our 
analysis. The recent LHC limit on gluino mass is not considered here. We will discuss it in 
detail in the next section. 

The results of our numerical analysis are displayed in Table (jl]). We study three different 



cases: 



1. In case I, we do not impose constraint OH]) and minimize R for fixed miQ = 10 TeV. The 
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Parameter 

77116 

my 2 
A /m 16 
m 10 /m w 

tan/3 
Am 2 H /m 2 w 



Case I 
10000 
34.05 
-2.29 

1.08 
51.39 

0.25 



Case II 
10000 

43.0 
-2.26 

1.11 
51.62 

0.28 



Case III 
15000 
51.07 
-2.29 
1.08 
51.39 
0.25 



Case IV 
20000 
62.65 
-2.45 
0.94 
54.81 
0.39 



R 



1.01 



1.04 



1.03 



1.02 



m?, 



X-1,2 

m-o 

A3, 4 



345.0 
49.6, 118.9 
6658.5, 6658.6 
119.6, 6650.6 



367.9 
53.2, 127.1 
6460.2, 6460.3 
126.5, 6452.6 



487.7 
73.6, 177.9 
9966.9, 9966.9 
183.7, 9953.2 



634.9 
101.4, 241.7 
17102, 17103 
242.9, 17097 



my 



in 



e L,R 



111 



*1,2 



mi. 



rrii 



9984.3, 9891.6 
9984.7, 10043 
9926.0, 7374.4 
9924.7, 10134 
2554.9, 3181.6 
2997.6, 3341.7 
3904.9, 7365.2 



9988.4, 9874.7 

9988.8, 10052 
9912.4, 7274.3 
9911.2, 10159 
2506.2, 3110.1 
2902.4, 3226.2 

3664.9, 7278.2 



14988, 14850 

14988, 15075 

14893, 11193 

14890, 15199 

3957.5, 4859.7 

4724.0, 5097.9 

6329.0, 11192 



20010, 19782 

20010, 20126 

19841, 14636 

19838, 20312 

5915.9, 6807.6 

6538.7, 6984.1 

7549.7, 14617 



m A 



126.0 
3463.6 
3441.4 
3465.5 



125.8 
4353.4 
4325.4 
4354.9 



127.8 
5664.2 
5627.4 
5665.3 



125.0 
8683.1 
8626.5 
8683.8 



3.09 x 10" 4 3.07 x 10" 4 3.06 x 10" 4 3.06 x 10" 4 

0.79 x 10" 4 0.79 x 10" 4 0.79 x 10" 4 0.79 x 10~ 4 

4.55 x 10" 9 4.22 x 10" 9 4.15 x 10" 9 3.96 x 10~ 9 



arf) 



BF(6 
BF(B 
BF(B S -> n+/f 



0.024 x 10" 10 
2737 



0.026 x 10" 10 
719 



0.008 x 10~ 10 
39866 



0.008 x lO" 11 
14115 



TABLE I. The benchmark solutions obtained for t — b — r Yukawa unification in 5'O(10)+HS model 
for positive \i. Different columns correspond to different cases discussed in the text. All masses 
are in GeV units. 



best unification found corresponding to R = 1.01. The solution predicts BF(B S — > 
/i + yU~) = 4.55 x 10 -9 which is slightly above the upper bound (jHJ). 



2. The dominant contribution to B s — > yU + /i~ in MSSM is proportional to m^ 4 27J where 
niA is the mass of pseudo-scalar Higgs. Using this fact, we simultaneously maximize 
rriA and minimize R in case II for = 10 TeV. As a result, we get a lower value of 
BF(£> S — > yU + /i~) which is consistent with experimental limit (JSJ). However one gets a 



S 



slight declination in Yukawa unification in this case. 



3. Case III and IV correspond to an obvious way of decreasing SUSY contribution to 
flavor violation namely, uplifting the SUSY scale. We take = 15 and 20 TeV which 
increase the masses of all SUSY spectrum including m^. The unification achieved in 
these cases is 97-98%. 

It is clear from the results of our analysis that a very good t — b — r Yukawa unification 
can still be achieved with "low" mie i n 5O(10)+HS model without violating the present 
experimental constraint on B s — > /i + fi~ . The calculated values of BF(b — > 57), BF(B — > tu t ) 
and Aa M shown in Table ([!]) are almost similar to their standard model values. The SUSY 
contributions to these processes are negligible due to the heavy sparticle spectrum one 
typically gets in t — b — r unified solutions for fi > 0. Among the other well known features 
of t — b — t unification with positive /i are relatively heavier Higgs ~ 125 — 130 GeV 
arising due to large m 16 and the condition A ~ — 2m 16 0, Q and pure bino like lightest 



neutralino which leads to the over abundance of the neutralino dark matter 11 . One would 



need additional mechanism, e.g. tiny R parity violation 29] to reduce this abundance. 



III. t-b-r UNIFICATION AND HEAVIER GLUINO 

As noted earlier, the t — b — r unified solutions for positive /1 generally require very light 
gluino mass < 500 GeV. The direct SUSY searches at the LHC has now excluded m 5 < 620 
GeV in S'O(10)+HS model [17|. As a result, the solutions displayed in first three columns 
in Table (JI]) are ruled out. It is recently shown that consistent t — b — r unification with 



heavier gluino can be obtained by increasing uiiq [18(. In fact one needs > 20 TeV to 
evade the present LHC bound on gluino mass e.g. case IV in Table ([!]). We propose here 
an alternate way to obtain heavier gluino in t — b — r unified solution without increasing 
mi6- As mentioned in Section [Til the ratio \Ay^\/\Ayf \ in Eq. ([7]) can get an additional 
suppression if m 4 ~ 2 < . This allows t — b — r unification with heavier gluinos. The required 
mass hierarchy m f ~ 2 < m~ b2 can be obtained if the appropriate squark splitting is introduced 
at the GUT scale. For example, consider SU(5) invariant boundary conditions: 

m % = m ! = m l = m 16 > 

m| = m 2 jj = m\ & + Am| . (11) 

The origin of such splitting in an 5*0(10) model will be discussed in the next section. As we 
will show later in Eq. fl3T| 132]) . Am| is allowed to take any values greater than — mf 6 . The 
choice Am| > raises the mass of one eigenstate of sbottom squarks, i.e. m~ b2 , compared to 
its value obtained with universal squark masses at the GUT scale. This leads to m~ t2 < m~ b2 
at weak scale which is the case of our interest. 

To quantify the effect of squark splitting, we study the above case through detailed 
numerical analysis. Fixing mi6 = 10 TeV we perform a random scan over all the remaining 
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soft parameters and tan/3. The analysis is performed for two scenarios: (1) only Higgs 

and (2) Higgs splitting + Squark splitting (HS+SS) with 
We employ the same numerical technique and apply all the theoretical 
and experimental constraints discussed in Section [TT1 The results of numerical analysis are 
shown in Fig. 03]). 



splitting, i.e. for Am| 
Am 2 H ,Am 2 s > 
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FIG. 3. Solutions of t — b — r Yukawa unification as a function of gluino mass for niiQ = 10 TeV. 
The figure in the left (right) panel shows the solutions obtained with only HS (HS+SS) at the GUT 
scale in SO (10). The vertical line corresponds to the lower bound on gluino mass in 5O(10)+HS 
model derived from the recent ATLAS data [13]. All the points shown are consistent with various 



phenomenological constraints discussed in the text. 



As can be seen from Fig. ([3]), the lower bound on gluino mass in SO( 10)+HS model 
rules out the Yukawa unification R < 1.08. In the presence of additional squark splitting, 
the unification up to 99% can be achieved without violating the present LHC limit on 
gluino mass. Also, a relatively heavier gluino up to 1.5 TeV is obtained assuming at most 
10% deviation in Yukawa unification without uplifting m 16 . It can be also seen that the 
number of valid solutions obtained with HS+SS are more compared to those obtained with 
only HS. In Fig. (@J, we show the ratio m^/m^ in HS and HS+SS models. The ratio 
substantially decreases in HS+SS model compared to its value without squark splitting. As 
we mentioned earlier in this section, this allows a heavier gluino in the spectrum without 
uplifting the SUSY braking scale. Note that with m 16 = 10 TeV, rrig > 1.5 TeV cannot 
be obtained for R < 1.1 even in the HS+SS model. This range of gluino mass is still in 
the reach of LHC and its future operations at = 14 TeV can significantly constraint the 
parameter space of HS+SS model if not rule it out completely. 



IV. t-b-r UNIFICATION AND REALISTIC 5O(10) 



We consider two categories of 5*0(10) models. One in which Higgs sector is extended 
to obtain realistic fermion masses and mixing and the other in which one introduces also 
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FIG. 4. Squark mass ratio m^/m^ obtained in HS (blue points) and HS+SS (red points) models. 



additional matter multiplet at Mqut- The former class of models lead to just so HS and 
the latter also to an independent squark splitting. We discuss them in turn. 



A. Just so HS in realistic 5O(10) 

5O(10) model containing a 10-plet of Higgs field 10h allows the following term in the 
superpotential: 

F 10 16 F 16 F 10 H , (12) 

where 16 p refers to the matter multiplet and Y w to the Yukawa coupling matrix in the 
generation space. The t — b — r unification follows under two assumptions: 

(A) Third generation of fermions obtain their masses only from Eq. f 1 1 2 p . 

(B) The MSSM fields H u and H d reside solely in 10 H . 

One however needs to violate both these assumptions in order to obtain correct masses for all 
fermions and non-trivial mixing among them. We estimate the effects of these violations on 
the t—b—T unification. We will take as an example a popular minimal renormalizable 5O(10) 



model 30] which is found to explain fermion masses and mixing in a number of situations, 



for instance, see 3JJ and references therein. The model contains a 126# field to generate 
neutrino mass and a 126# to preserve supersymmetry at the GUT scale. In addition, it has 
a Higgs transforming as 210# representation of 5O(10) which breaks 5O(10) to MSSM. 

Due to the presence of additional Higgs fields particularly 210 H , SU{2) L doublets resid- 
ing in various Higgs fields mix with each other. This mixing plays an important role in 

But as we show below, this 



30 



generating the right type of the second generation masses 
mixing also generates just so Higgs splitting required to obtain t — b — r unification if the soft 
masses of the 10#, 126// and 210# fields are non-universal. On the negative side, the pres- 
ence of 126 h and the Higgs mixing also lead to departures from the exact t — b — r unification. 
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The 10h, 126#, 210# fields each contain four down-type and four up-type standard 

model doublets. They mix and produce four mass eigenstates denoted as h^ ,d : 



K 4 = U a f4> U a' d , (13) 

where <^ d are components of bi-doublets in 10#, 126#, 126#, 210#. One assumes that 
through fine tuning only H u = h\ and H d = hf remain light. Consider now the soft mass 
terms of various SO(10) Higgs prior to the 5*0(10) breaking: 

Vsoft 3 m 2 w 16^16 F + m 2 w 10^ H 10 H + m? 26 (126^126 H + 126^126^) + m 2 10 210^210^ . (14) 



Here we have assumed equal masses for the 126# and 126# fields to avoid non-zero D-term. 
The masses of the other Higgs multiplets are taken non-universal. Substitution of Eq. ( |T3|) 
in Eq. (HU) leads to 



V soft 3 m 2 l6 {&Q + rfU + D j D + D~L + E*E) 

+ Ko|Or i r + mL 6 (|0^| 2 + |0 3 ^| 2 )+m2 10 |0 4 u 1 | 2 )^ 

+ (u -> d). (15) 
The above equation leads to the following boundary conditions at the GUT scale 

2 2 2 2 2 2 

wiq = rrijj = = = = m 16 , 
m H u , d = m ?26 + \On d \ 2 (m 2 w - m 2 26 ) + \0^ d \ 2 (m 2 2W - m 2 26 ) . (16) 

It is seen that the Higgs mixing generated through Eq. ( {TBI) has produced the desirable 
splitting only among H ud masses without splitting squarks from each other unlike in case 
of the popular D-term splitting. 

Let us now look at the impact of Higgs mixing on the t — b — r unification. The presence 
of m> H field modifies Eq. fllj) to 



r 10 16 F 16 F 10^ + F 126 16 F 16 F 126^, (17) 



where Y 126 additional Yukawa coupling matrix. By substituting Eq. (113]) in Eq. (Tr7|) one 
arrives at the charged fermion mass matrices: 

M d = v d (Y w O d 1 + Y 126 O d 1 ), (18) 
M l = v d (Y l0 O d ll -3Y lw O d 2l ), (19) 
M u = v u {Y w O u 11 + Y l2 ,O u 21 ) 1 (20) 

where v U)d denote the vacuum expectation values of the neutral component of H Utd . We can 
go to a basis with Y w diagonal. Neglecting the contribution of 126# for the time being, one 
has 

Qd 

Vb = Vr = y3vt- (21) 
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Thus one source of the departure from the t — b — r unification is the ratio O^/O^. The 



parametric form of the matrices O u,d is worked out 32|, [33] in the model under consideration 



and we closely follow the notation in [33|. This is based on the Higgs superpotential 

W H = M w 10 2 H + M 210 210 2 H + M 126 1 26 H T26 H 

+ A 210^ + r) 210 H 126 H 126 H + 210^10 H (al26^ + aY26 H ). (22) 

The Higgs mass matrices and hence O u,d follow from the above superpotential after the 
50(10) breaking and in the most general situation with 50(10) breaking to standard model 



one obtains (see, Eqs. (C18, C19) in [331 ]) 

H d = N d (-^-<f>i - V6-(2x -l)(x + l)^<p d 2 - V6-(3x - l)(x 3 + 5x - + a—p' 3 <p d 
\x—l r) P3 V rn^ 

H u = N u ( -^-0? - VE-{2x - l)(x + 1)-^ - V6-(3x - l)(x 3 + 5x - 1)0^ - a— y'^l 
\x - 1 T] P3 V rn<t> 

(23) 

where N u ^ are overall normalization constants, x is an arbitrary parameter andp3,p 3 ,]?5 are 
polynomial in x. The expressions of elements of O u,d can be read from the above equation. 
The ratio 0^/0^ which measures deviation from t — b — r unification can be exactly or 
nearly one in a number of situations. Obvious case is the limit N u = Nd « 1 corresponding 
to the situation in which extra Higgs fields' 02 3 4 contribution to H ud are sub-dominant. 
In general, various 0^ and 0f are components of SU(2) L x SU(2) R bi-doublets residing in 
10h, 126#, 126 h, 210h representations. They are thus distinguished by the SU(2)r group. 
One therefore automatically has O u = O d as long as SU{2)r is unbroken. This happens 
331 ] for x — 0. One can then show from Eq. (T23]) that O u = O d in this case. Another 
interesting limit corresponds to choosing a = a in Eq. (|23D. In this limit, O u ^ O d but 
still 0" x = Of x thus one obtains exact y& = yt in Eq. (121!) even in the situations where 
SU{2) R is broken. This case also corresponds to = —0\ x and thus HS also vanish in 
this limit as follows from Eq. (TT6]) . But mild deviation from the limit a = a can generate 
sizable HS and approximate t — b — r unification for large ranges in other parameter. This is 
illustrated in Fig. (jSD where we show ^ and l^ 3 -] as a function of x for two specific values 

P*d ^41 

of e = l/2(a — a) = 0.1, 0.2. The remaining unknown parameters appearing in Eq. f|23|) 
are equated to 1. It is seen that for most values of the unknown parameter x one obtains 



34 



as given in Eq. (fib]) . 



almost exact t — b — r unification, i.e. — 1 and non-zero HS 

Another threat to the t — b — r unification comes due to the presence of the Yi 2 q Yukawa 
couplings in Eq. (1181) . This effect is somewhat model-dependent and we estimate it by 
specializing to the case of the second and third generations. We can write the charged 
fermion mass matrices as 

%Jf ( h 2 + x 2 x \ ( h 2 -3x 2 -3x 

M d = v d \ ; Mi = v d \ 

\ x h 3 + x 3 I \ —6x n 3 - 3x 3 

Mu = Vu 91l ( h + 3X2 SX V (24) 
u Of 1 \ sx h 3 + sx 3 J ' 
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FIG. 5. The Higgs mixing parameters jf^ (continuous lines) and \j4r\ (dashed lines) as a function 

a 41 

of x in the minimal SUSY 5*0(10) model. The blue (red) line corresponds to parameter e = 0.1 
(0.2). 



QU Qd 

where s = 7^-7^-- Here /i 2 ,3 refer to elements of the diagonal YiqO^ and x 2 , x, x 3 to that of 
21 O u 

symmetric O^Y^- Several of these can be determined from the known masses and mixing. 
Approximate t — b — r unification is obtained with the hierarchy x, x 2 ,x 3 *C ^3- Assuming 
h 2 <C X2 then leads to the desirable mass relation 3m s = m^. Given this hierarchy one finds 
for real parameters 

m s T/ M % m b m c 

x 2 ~ — ; x ~ m b V cb {l - s); s~ 

mft m s mt 

The ratio x 2 /x 3 remains undermined. If type II seesaw dominates and Yyi§ is to provide the 
neutrino masses and mixing then x 2 ~ £3 [351 ] . In this case, one finds 



2/6 



^3 (1 + ^) (i + o(k 2 6 )), 
(i + o(v; 2 b )) 



On V m < 



where we have used Eq. (pQ). Using value for the mass ratios at the GUT scale for tan (3 = 50 



36 



^ rs 0.016, ^ ~ 0.0023 already implies about 6% departure from 6 — r unification. 
The m s ,m c dependent terms in Eq. (]2"5"|) get suppressed if one assumes x 3 <C x 2 . In this 
case, one can still obtain b — r unification and reproduce the second generation masses and 
Cabibbo mixing. The Higgs mixing factor can be nearly one as argued before. It is thus 
quite plausible that one can obtain almost exact Yukawa unification not just in simplified 
but also in more realistic GUT based on 5*0(10). 
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B. Squark splitting and t — b — r unification 



The squark splitting can be induced by adding new matter fields having the same quan- 
tum numbers as some of the squarks and letting them mix with the normal squarks. The 
minimal possibility at the 5*0(10) level is introduction of a 10m field. This contains fields 
transforming under an SU(5) subgroup of SO (10) as 5' M + 5 M . Of these, the 5 M can mix 
with 5a/ contained in the matter multiplet 16m of 50(10). This mixing can generate the 
squark splitting. Let us discuss details within a specific model which has been studied ex- 



tensively |37|, |38[ for several different reasons. The model contains three generations of 16m 
and three copies of 10m- We will explicitly consider only the third generation and a 10m in 
the following. The Higgs sector consists of the usual 10# supplemented by 16# + 16h and 
45// + 54#. The 16// is introduced to break the B — L gauge symmetry and 45// + 54// are 
needed to complete the breakdown of 5O(10) to SM. This model serves as a good example 
in which (a) independent squark and Higgs splitting can be generated and (b) there exit 
ranges of parameters for which t — b — r unification is approximately maintained. We shall 
discuss only a part of the superpotential relevant to describe Higgs and squark mixing, see 



381 ] for a general discussion of the model. The superpotential of the model can be divided 
in two parts one describing matter-Higgs interaction and the other describing Higgs-Higgs 
interactions: 

W M = F16m16m10// + F16m10m16// + M10 m 10m, 

W H = M 16 T6//16// + M 10 10//10// + H16 H 16 H 10 H + H'TE H 1Q H 10 H . (26) 

The above superpotential is designed to respect the matter parity under which all the matter 
(Higgs) fields are odd (even). This symmetry is essential for preventing renormalizable 
baryon and lepton number violating terms. Scalar components of none of the matter fields 
acquire vacuum expectation value (vev) and thus matter parity remains unbroken. Only 
fields appearing in the above superpotential and acquiring the GUT scale vev are 1// + 1// 
contained in 16// + 16// of 5O(10). Thus after the GUT scale breaking, above superpotial 
maintains invariance under the SU(5) subgroup of 5O(10). As a result, the mixing between 
the following SU(5) components is allowed and arise from Eq. ( J26l) : 



5 M \ _ r>( Q \ { 5 m \ . / 5# \ _ , . / 5// \ / 5^ 





where 




m ;r , (27) 



cos j -sinj 
sin j cos j 



The fields with (without) prime are component of the original 10 (16 + 16) of 5O(10). It is 
assumed that fields labeled with superscript I are kept light by fine tuning and those with 
the superscript h pick up masses at the GUT scale. The mixing angles appearing above are 
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related to parameters in Eq. (|26|) and are explicitly given as 

Fv 16 



tan# 
tan 27 



M ' 

2v 16 (HM w + H'M 16 ) 
(M? 6 -M? Q )+vl 6 (H*-H'Z) 
2v 16 (H'M 10 + HM w ) 



tan2 ^(M^-M l2o )-^-^)' ™ 

where vi e is a vev of SU(5) singlets residing in 16# + 16#. The light fields transform under 
S77(5) as 5h + 5# of Higgs and 10m + 1a/ + 5a/ which together makes 16a/ and 10# of 
SO (10). The effective Yukawa couplings of these fields can be worked out from Eqs. (12"B"| 
|27|) in a straightforward way. One finds 

W eff = Ys 5 10 M 10 M 5 l H + (Ys 7 c e + Fs e c^lO M 5 l M 5 l H , (29) 

where Sj = sin j and Cj = cosj. 

The Higgs mixing as given in Eq. ( 1271) leads to both the squark and Higgs splitting 
through SO(10) invariant non-universal boundary conditions. Consider the following soft 
terms: 

V sof t 3 m? 6 16 M 16A/ + m? 10 M 10 M + ro? 6 16^16* + m%ltf H lQ H 
= rr4(10 M 10 M + 51,5a/ + 1 a/ 1a/) + m 2 10 (5%5' M + 5 M 5 M ) 
+ m? 6 (10\jl0 H + %5 H + l\jl H ) + mf (5j5^ + 5^). (30) 

Substitution of Eq. (1271) in Eq. (1301) leads to the following boundary conditions for the soft 
mass parameters of squarks and Higgs: 

m \ = m }j = m \ = m i6' 

2 2 2,2/2 2 \ 

2 /2 2 1 12 2 

m // d = m l6 C 7 + m lO S 7' 

2 /2 2 1 /2 2 /ni\ 

m Z/„ = m 16 C 5 + m lQ S 5 ( 31 ) 

resulting in the following squarks and Higgs splitting 

A 2 2 2 2/ 2 2 \ 

Am 5 = - = s e (m 10 - m 16 ), 

Am^ = - m 2 Hu = (s* - s])(m% - mf 6 ). (32) 

Thus model under consideration simultaneously generates independent mixing among 
squarks and Higgses that result into squark splitting and Higgs splitting respectively. The 
mixing angles which generate these splitting also lead to departure from the exact t — b — r 
unification as before. But the exact b-r unification follows for arbitrary values of these 
mixing angles. Even t — b — r unification also holds approximately in limiting cases, e.g. 
Fse <C 1 and H = H'. Another interesting limit corresponds to ss ~ s^cg and Fsg <C 1. In 
this limit one gets yj, = y T ~ yt and simultaneously non-zero splittings, Am 2 H pa s 2 Am 2 s for 
m w = m' w and m 16 = m' 16 in Eq. (1321 . This relation automatically implies Am^ > as 
required for the REWSB when Am| is chosen positive to raise the gluino mass limit. 
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V. SUMMARY 



We have addressed two important issues in this paper in the context of the t — b — r 
unification in 50(10) broken to MSSM with a positive \i parameter. The existing analysis 
[9 -12 1 have either assumed only HS with degenerate squarks at the GUT scale or a D-term 



splitting in which case squark splitting is correlated to the HS. This scenario appears to be 
ruled out save for very high SUSY scale m 16 around 20 TeV [l8[ . Detailed phenomenological 
analysis presented here shows that independent and positive squark splitting Am|, Eq. (fil~j) 
can change the allowed gluino mass and scenario can be compatible with t — b — r unifica- 
tion R ~ 1.01 and the recent ATLAS bound m g > 620 GeV. Moreover, R < 1.1 requires 
rrig < 1.5 TeV. Thus viability of the t — b — r unification can be tested in future at LHC with 
y^s = 14 TeV. The other issue addressed here concerns the origin of just so HS. While just 
so HS is introduced as a phenomenological parameter in many works |9j-[l2j, its origin is not 
justified in most of the existing analysis, see however 00. We have taken here a concrete 



realistic model 30] used to understand fermion mass spectrum and shown within it that 
just so HS is an automatic consequence of the non- universal boundary conditions and Higgs 
mixing. It is also shown that one can obtain sizable just so HS and retain almost exact 
t — b — t unification in this realistic model. Independent squark splitting required to relax 
the gluino mass bound is also shown to follow in an extended model in which squarks mix 
with additional matter multiplet. One may conclude based on the analysis presented here 
that t — b — r unification with positive \i is still phenomenologically viable and theoretically 
well-founded. 
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